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Abstract. Influence of the damping on collective oscillations of a one-dimensional 
trapped Bose gas in the mean field regime has been studied. Using the 
phenomcnological damping approach developed by L.P. Pitaevskii, modified variational 
equations for the parameters of the condensate wave function is derived. Analytical 
expressions for the condensate parameters in equilibrium state have been obtained. 
Bistability in nonlinear oscillations of the condensate under periodic variations of the 
trap potential is predicted. The predictions of the modified variational approach are 
confirmed by full numerical simulations of the ID GP equation with the damping. 
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1. Introduction 

The dynamics of a one-dimensional trapped ultra-cold Bose gas has attracted 
considerable attention for last years pp. Recently has been achieved regimes where ID 
Bose gas with a zero-range two-body interaction in the Tonks -Girardeau (TG)regime 
El H] with fermionic behavior becomes visible. An investigation of the transition 
from the mean field regime to the TG regime was performed in [5]. Experimentally ID 
regime has been realized in [Hj. 

Measurements of the collective oscillations of such a system should give a lot of 
information about the BEC dynamics. In particular this is important in the analysis of 
the condensate dynamics in a magnetic waveguide, being a fundamental atom optical 
element [7j. Performed by this time theoretical descriptions have mainly dealt with 
conservative systems, e.g. see [Hj, where collective oscillations of a quasi ID Bose- 
Einstein condensate (BEC) in the low and high density regimes were investigated. The 
damping of the radial BEC oscillations in a cylindric trap, connected with the parametric 
resonance and leading to the energy transfer from collective oscillations to longitudinal 
sound waves has been studied in 0). The dissipation inheres in real systems. So it is 
of interest to investigate theoretically effect of damping on collective oscillations of a 
one-dimensional trapped repulsive Bose gas. 

We consider here the problem using the phenomenological approach developed by 
L.P. Pitaevskii jTU] and employed later in [TTJ 

2. The model 

The dynamics of a trapped one dimensional repulsive Bose gas with the damping is 
described in the framework of the ID Gross-Pitaevskii equation fHJ ITT] 



with the total number of atoms N = J \(p\ 2 dx. This equation is obtained in the case 
of a highly anisotropic external potential under the assumption that the transversal 
trapping potential is harmonic: V(y,z) = muj_(y 2 + z 2 )/2 and u>± ^> uj x . Under such 
conditions we can consider the solution of 3D equation to have the form U(x, y, z; t) = 
R(y, z)(f)(x, t) where = mwi exp(—mu}j_p 2 /h)/(nh). Averaging in the radial direction 
(i.e. integrating over the transversal variables) we come to equation describing the 
dynamics of the gas in longitudinal direction. The potential V(x, t) is assumed to be 
V(x, t) = mui 2 x 2 F{t) , where F(t) describes the time dependence of the potential, 
which we consider here as F(t) = 1 + hsin(wt). The effective one dimensional mean 
field nonlinearity coefficient g^o = 2ha s u±, where a s is the atomic scattering length. 
a s > corresponds to the Bose gas with a repulsive interaction between atoms and 
a s < to an attractive interaction. 




(1) 
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We have the following dimensionless form of the equation (JTJ 
1 x 2 

= i7(-^*x + + ~ /#) = Rty, r), (2) 



by setting: t = u x t, I = Jh/(mw x ), x = x/l, ip = J 2\a s \uj_ L /uj x (f) } with g = ±1 for the 
repulsive and attractive two-body interactions respectively. 

3. Variational analysis 

To describe collective oscillations of a Bose gas under damping we employ the variational 
approach. For the wavefunction ip(x,t) we use the gaussian trial function 

x 2 ih(t)x 2 

4>(x, t) = A(t)exp(-^-^ - - i<p(t)), (3) 

where A, a, b and ip are the amplitude, width, chirp and linear phase, respectively. 
Equation (J2J can be obtained from the variational equations 
dL d dL d dL 5L R n 

I = (4) 

dip* dxdip* dtdipt 5ip* K ' 

where L is the Lagrangian density, L = L(x,t), of a conservative system, given by 

l = - m - \m 2 - (y^w - »m 2 -\w (5) 

and L R is defined as dL R /dip* = —R(ip, if)*), where R(ip, ip*) is the right side of equation 

Inserting trial function^ into equation © and averaging it as 

L(x,t)dx (6) 

we obtain the averaged Lagrangian of the conservative system in terms of the trial 
function parameters: 

L A 2 a\ A2 A 2 A 2 a 3 b 2 

= = + A 2 a<p t 

7T 4 ^ 4a 4 

FA 2 a 3 gA 4 a a2 

1 \xA 2 a. (7) 



4 2^2 

Using equation (HJ) and its conjugate, we obtain a system of equations for the 
variational parameters rji [TBI H4j: 

dL d9L f dx(R d -f + R^). (8) 



dr]i dtdrjit J % % 
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Inserting (j3J) and (J7J) into equation (jSJ) we derive the following system of ordinary 
differential equations (ODE): 

d(A 2 a) _ 7A 2 1 A 2 o?b 2 1 FA 2 a 3 
dt ~ ~2a~ + 2 + 2 + 

+V2^ygA 4 a - 2<y/j,A 2 a, 

d(A 2 a 3 ) a2 o 7 A 2 a 3 7 A 2 a 5 6 2 
v ; - -2A 2 a 3 b - -— + — — + 



dt 2 2 

iA 4 a 

72 



H 1 7= 27/iA a 



db 2 7 6 1 L2 ^ c/^ 2 
-77 = ^ r + & 2 + F 



dt a 2 a 4 v/2a 2 ' 
<fy> lb 1 5gA 2 

^ = -y + 2^ + vf-^ (9) 

This system of equations can also be obtained from the modified form of the 
conservation law for the number of atoms and by using the moments method, as shown 
in Appendix. Let us rewrite the system using the notations x = a 2 ,y = A 2 : 

x t = —2x6 — 7 + 7X 2 6 2 + j Fx 2 



V2 



Vt = yb + — + —7=- - 2<yny, 
x 2V2 



6, = ^_J_ +62 + i ,_ jy. (10) 

x x 1 V2x 
This ODE system is the main result of this section. 

4. Numerical simulations 

We have carried out a series of time dependent simulations of the system within 
the variational approach using equation (JIUj) and also by performing exact numerical 
calculations using equation (J2J). In our numerical calculations we discretize the problem 
in a standard way, with the time step dt, and spatial step dx, so ipj approximates 
ip(jdx,kdt). More specifically we approximate the governing equation (J2J) with the 
following semi-implicit Crank-Nickolson scheme using split-step method [12J 

*(^ +1 -^ +§ ) *(^ +§ -^) 
dt dt 

— dt = 2 ^ j ^ j ' 



where H x = 0.5(1 + i 7 )(x 2 F(t)/2 + g\ij(x,t)\ 2 - //), H 2 = -(1 + ij)d 2 /dx 2 , H 3 = H x 
and ip k+ 3, ip k+ 3 are defined so that 

^£izM = \h^ + ^), (12) 
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(13) 



k+l 




iF 3 (^ +1 +^ +l )- 



(14) 



Solving ()12|) and (|13p first at time tk = kdt we produce intermediate solutions ip k+ z and 
?/> fc+ t. The final solution for one time step dt is obtained from (|14jl . 

The results of numerical simulations of both PDE and ODE models are presented 
below. 

Figure ^ shows time evolution of the width and the norm of the condensate for the 
case of F(t) = 1, i.e. when the external potential has no time perturbations. The figure 
is presented to compare the character of the width oscillations and the behavior of the 
norm depending on the values of the chemical potential /J and the dissipative constant 
7. As seen, the frequency of the oscillations does not depend on the dissipative constant 
at all and depends weakly on the chemical potential, e.g. when fi — 2, w = 1.808 while 
when fj, = 3, wo = 1.77. For greater values of 7, the oscillations damp faster. We see 
that the ODE leads to an equilibrium state the norm of which is 3-4 percent less than 
that of the PDE results. 

From figure [T] it can be seen that damping process of the condensate eventually 
leads to the equilibrium state. This equilibrium state can also be obtained analytically 
by solving equation pOjl . Taking into consideration that in an equilibrium state xt = 0, 
y t — and 6 = the following expressions can be obtained: 



Putting the value of chemical potential e.g. for /x = 2 we find a = 1.383, A = 1.403 and 
N = A 2 ay/ir = 4.825 which are confirmed by the PDE results. 

In figure El we can observe an interesting behavior of the norm. If an external 
perturbation is applied to a trapped BEC which is already in the equilibrium state then 
the norm of the condensate starts decreasing. In the figure the frequency of the periodical 
trap perturbation is taken to be equal to the eigenfrequency of the system determined 
from the figure ^ and the amplitude of time perturbations is taken as h — 0.06. We 
see that for smaller values of the dissipative constant 7 the norm goes farther from the 
equilibrium state. That is, the dissipative constant is not simply the quantity which is 
responsible for diminishing the norm, but it is the constant keeping a condensate in the 
equilibrium state. 

The width dynamics under main resonance with the initial wave packet taken in 
the equilibrium state is depicted in figure El As shown, in contrast to the norm the 
width oscillates near the previous point. As seen from the figure the amplitude of the 
width oscillations becomes stable by the time t = 140. 



A 2 



Afi + y/TEj? + 60 

10 

4^2 2^2 



(15) 



5 5a 2 
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Performing ODE and PDE calculations for the frequencies which lie around the 
eigenfrequency of the BEC and measuring the amplitudes of the above-said stable 
oscillations we have plotted the values of the oscillation amplitude as a function of 
the frequency of the periodical trap perturbations with h = 0.03 and h = 0.06 in the 
cases 7 = 0.01 and 7 = 0.005 in figure |U 

We see that e.g. when h = 0.06 and 7 = 0.005 the highest amplitude of oscillations 
is driven in the trap perturbation with the frequency w = 1.89 which is more than 
the eigenfrequency wq = 1.806. Bistability appears with the smaller values of 7 in 
the vicinity of this critical frequency. For w = 1.88 we observe large oscillations with 
a osc = 1.6, while at w = 1.91 we observe much smaller width oscillations with a osc = 0.6. 
It can be seen that with the growth of trap perturbations the value of the critical 
frequency becomes greater. 

Let us estimate parameters for the experiment. The magnetic trap can be taken 
with parameters u>j_ = 2tt x 40 Hz, uo x = 2tt Hz, and the number of atoms of 85 Rb 
iV = 0.23 x 10 4 . For the external field B = 161.57, a s = 0.3nm (repulsive gas). Then, 
the eigenfrequency of the harmonically trapped BEC 27rz/ = 1.808 x uj x = 11, 36 Hz. By 
applying the external perturbation F(t) = 1 + 0.06 sin(27rz/t) with 2-nv = 11.36 Hz to the 
trapped BEC one will observe decreasing of the number of atoms by 20, 13 and 7 percent 
for 7 = 0.1, 0.2 and 0.3 respectively. For the external perturbation with 2ttu = 11.8 Hz, 
large oscillations will be observed with a osc = 1.6/ = \.Q\jTi/ {muj x ) = 17.4/im, while at 
2ttv = 12 Hz one will observe much smaller oscillations with a osc = 0.6/ = 6.5/im. 

5. Conclusions 

In this paper we study the collective oscillations of a quasi one dimensional Bose gas 
in the presence of dissipative effects. The modified Gross-Pitaevskii equation in the 
framework of the phenomenological approach [TUl ITT] has been used. To describe 
evolution of oscillations we employ the modified variational approach taking into account 
the dissipation. We confirm results for the oscillations damping obtained from the 
system of equations for the wave function parameters, by the moments method and 
direct numerical simulations of the full GP equation. 

The calculations show that damping oscillations of a BEC which eventually comes 
to an equilibrium state occur with the eigenfrequency which does not depend on the 
value of the dissipative constant 7 and only depends on the chemical potential /1. 

The expressions for the width and the norm of a condensate in an equilibrium state 
have been derived analytically. 

We study the main resonance in condensate oscillations. For a BEC in an 
equilibrium state it is shown that making the external potential oscillate leads to 
decreasing of the condensate norm and the norm begins oscillating around the point 
which is less than the previous stationary one in the equilibrium state, whereas the 
condensate width oscillates near the previous point. 

We have also shown that in resonances the bistability appears with smaller values 
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of the dissipative constant in the vicinity of the critical frequency which is above the 
eigenfrequency of the BEC. 
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Appendix 

Differentiating the norm of the condensate by time we can obtain the first equation of 
the system (jTUJl as 

" ' \^\ 2 dx = l(^ + ^t)dx. (A.l) 



dt dt 

Inserting ip t from equation (J2J) into (|A.1|) we derive the modified form of the 
conservation law for the norm of the condensate [29] . 
dN 



7 / M dx +jF J x 2 \ip\ 2 dx + 
+2 7 # [ \ip\ 4 dx-2~ffi [\ip\ 2 dx. (A.2) 



Substituting the Gaussian trial function (jHJ) into this equation we obtain equation (fTUk ). 

Equations (fTUb) and (fTUb) can be obtained by the moments method. We obtain 
equation (fTUb ) by calculating 

' ''" ■ \i)\x 2 i) + i)*x 2 i) t )dx. (A.3) 



dt 

The substitution of the ipt from equation (J2J) leads to the following equation: 

d <c x 2 ^> r r 

— = i J x(ip* x ?p - ip x ip*)dx + 7j x 2 \ip x \ 2 dx - 

-7 J \ip\ 2 dx + jF J x A \ip\ 2 dx + 

+2 7 # j x 2 \ip\ 4 dx - 2 7 /i J x 2 \ip\ 2 dx. (A.4) 

Substituting trial function Q we get equation (fTUb ). Equation (fTUb ) is derived by 
calculating 

d -^dT^=-jt! r ^ dx= 

(ipt^xx + ip*iptxx)dx. (A. 5) 



Integrating by parts and substituting ip t from equation (j2J) we obtain 

— ( \ip x \ 2 dx = iF f x(ip*ip x — ipipDdx + 



dt 

+ig [^ 2 x -^f)dx + 
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+7 / \^ xx \ 2 dx-^F J \i)\ 2 dx + 

+ 1 F J x 2 \tfj x \ 2 dx + 4 ig J \tfj\ 2 \^ x \ 2 dx + 

+19 J {r 2 € + ip 2 r x 2 )dx - 2 7 /i J \^ x \ 2 dx. (A.6) 

Inserting the Gaussian ansatz into this equation and using results (fTTik ) and fTOb ) 
we obtain equation (flUb ). 
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Figure 1. Dynamics of the norm (a) and the width (b) of the repulsive BEC in a 
harmonic trap without perturbation. The upper and lower diagrams are plotted for 
the cases u = 3 and u = 2 respectively. Solid and dotted lines show PDE and ODE 
results. 



Figure 2. The behavior of the norm of the trapped BEC when the trap starts 
oscillating with the amplitude h = 0.06. Before the trap oscillations BEC was in 
the equilibrium state. We compare the PDE and ODE system simulations. The solid 
lines stand for full numerical simulations of the PDE, while the dotted lines represent 
the ODE results. 



Figure 3. BEC width versus time t in the main resonance in the case h = 0.06, \i = 2, 
7 = 0.01. Solid and dotted lines show the PDE and ODE results respectively. 



Figure 4. Oscillation amplitude as predicted by the PDE (scatter) and ODE (line) 
models. Here 7 = 0.005 (solid lines and squares) and 7 = 0.01 (dotted lines and 
circles). The upper two lines are for the case h = 0.06, the lower two are for the case 
h = 0.03. The chemical potential is equal to 2. 
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